Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



A multi-Frequency Study of The Radio Galaxy NGC326. 

I. The Data 

M. Murgia 1,2 , P. Parma 1 , H. R. de Ruiter 1 ' 4 , M. Bondi 1 , R. D. Ekers 5 , R. Fanti 1 ' 3 , and E. B. Fomalont 6 



Istituto di Radioastronomia del CNR, Via Gobetti 101, 1-40129, Bologna, Italy 
Dipartimento di Astronomia, Universita di Bologna, Via Ranzani 1, 1-40127 Bologna, Italy 
Dipartimento di Fisica, Universita di Bologna, Via B. Pichat 6/2, 1-40127 Bologna, Italy 
Osservatorio Astronomico di Bologna, Via Ranzani 1,1-40127 Bologna, Italy 
Australia Telescope National Facility, CSIRO, P.O. Box 76, Epping, NSW 2121, Australia 
National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 2293 



Received; Accepted 

Abstract. We present the results of a multi-frequency study of the inversion symmetric radio galaxy NGC326 based 
on Very Large Array observations at 1.4, 1.6, 4.8, 8.5 and 14.9 GHz. The morphological, spectral and polarization 
properties of this peculiar object are studied at different levels of spatial resolutions. The interpretation of the 
data will be discussed in forthcoming papers. 
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1. Introduction 

The radio source B2 0055+26 was identified with the 
elliptical galaxy NGC326 during a program directed at 
the identifications of the optical counterparts of radio 
sources selected from the B2 catalogue ( Colla et al. 1975] ). 
Basic properties of NGC326 are reported in Table 1. The 
radio source was first mapped by Fanti et al. (1977) at 
1.4 GHz with the Westerbork Synthesis Radio Telescope 
(WRST) in the snapshot mode. Because of the interesting 
structure seen, this observation was followed up by a 
more sensitive WRST observation at 5 GHz with an 
angular resolution of 6" x 13" QEkcrs et al. 1978|). At this 



called 'X'-shape, is characterized by two separate low 
brightness wings, almost perpendicular to the currently 
active lobes. There is a broad consensus concerning the 
dynamical interpretation of the 'C'-shape. This kind of 
distortion is cither caused by the translational motion of 
the galaxy through the intergalactic medium (wide and 
narrow angle tails) or by the orbital motion of the galaxy 
around a nearby companion ( Blandford fc Icke 1978 ). 
The explanation of the 'Z '-shape and 'X'-shape is 
still uncertain. Along with the secular jet precession 



resolution the radio source was found to be composed of 
two curving tails with a striking 180° symmetry. NGC326 
was recognized as one of the most spectacular examples of 
inversion symmetric radio galaxy. The authors explained 
the unusual structure as being due to the slow precession 
of the radio jets during the lifetime of the radio source. 
Rees (1978) suggested that the beam direction precesses 
due to a realignment caused by the accretion of gas with 
a different angular momentum direction with respect to 
the central black hole axis. 

Distorted radio galaxies can by classified according 
to two broad schemes: the mirror symmetric ('C'-shape) 
and the inversion symmetric ('Z'-shape). A sub-class, 
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(Rccs 1978), other two alternative scenarios have been 
proposed: the sudden realignment of the jet (Wirth et al. 
1982) and the buoyancy of the material from the lobes 
( |Worrall et al. 1995| ). 

Battistini et al. (1980) found that NGC326 is a double 
system composed of two nearly equally bright elliptical 
galaxies in a common envelope ("dumbbell" galaxy). 
Wirth et al. (1982) investigated the connection between 
the optical and the radio morphology for a sample of ^100 
radio emitting dumbbell galaxies. They found about a 
dozen of sources with 'Z' or 'X'-shape and suggested that 
the presence of another equally massive galaxy within 
10-30 kpc of a radio galaxy would strongly influence the 
jet properties. Particularly, they argue that in the case 
of bound circular orbits the continuous tidal interaction 
can produce mirror-symmetric wiggling jets as previously 
suggested by Blandford & Icke (1978). In the case of 
unbound orbits an impulsive interaction between the two 
galaxies can lead to an inversion-symmetric radio source 



2 



M. Murgia et al.: Multi- Frequency Study of The Radio Galaxy NGC326. I. The Data 



like NGC326. The 'Z'-shape occurs when the duration 
of the impulsive torque is longer than the jet outflow 
time to the lobes, while the 'X'-shape results when the 
crossing time is so fast, as compared to the jet outflow 
time, as to cause an instantaneous change in the jet 
direction. Wirth et al. (1982)| studied in detail the case 
of NGC326, using a new optical image. By examining 
the isophotes they found a brightness difference <1 mag 
between the two galaxies and very little deviation from 
pure ellipticity. Indeed they proposed NGC326 as a 
prototype for the class of undisturbed dumbbell galaxies. 
They also pointed the attention to the quite high velocity 
difference between the two galaxies: Av ~ 750 krns" 1 
(Sargent 1973). In the meantime the first 20-cm Very 
Large Ar ray (VLA) image of the ra dio source became 
available flFomalont 198l| , [Ekers 1982;). In the lig ht of the 
new optical and radio data, Wirth et al. (1982) reviewed 



Table 1. Basic properties of NGC326. 



the morphological classification and the interpretation of 
the source made by Ekers et al. (1978). They suggested 
that NGC326 describe an 'X'-shape in which the old jet 
direction was from NE to SW and the actual jet direction 
is from SE to NW. Moreover, both the old and the new 
lobes show a 'C '-symmetry but with a difference in the 
position angle of about 120° , suggesting that the radio 
galaxy has abruptly changed its velocity direction during 
its transient interaction with the other galaxy. In the 
Wirth et al. (1982) picture the velocity of the jet flow 
cannot be much grater than the Av of the galaxies. More 
recent redshift measures give Av = 509 ± 44 km s 



( Davoust fc Considere 1995| ) and Av = 549 ± 38 km s 
(Werner et al. 1999|). 



l 

-l 



NGC326 is the brightest member of a small group of 
galaxies, Zwicky 0056.9+2636. Werner et al. (1999) mea- 
sured the redshift for eight galaxies of the group (including 
NGC326) confirming the presence of a cluster with a mean 
redshift of z mean = 0.0477 ±0.0007 and line-of-sight veloc- 
ity dispersion a z = 599t^° krns" 1 . Werner et al. (1999) 
showed that the brightest of the two optical galaxies (the 
radio source) is a slowly-moving member of the cluster 
while its companion has a velocity of about 500 km s _1 
relative to the cluster velocity. 

The first detailed studies of the intergalactic medium 
surrounding NGC326 were done by Worrall et al. (1995) 
and Worrall & Birkinshaw (2000). The region containing 
the source was imaged in soft X-ray with the ROSAT 
PSPC. Surprisingly, they discovered that the galaxy is 
embedded in a bright asymmetric X-ray-emitting gaseous 
medium with properties more typical of a cluster, rather 
than of a group. They report a gas temperature of kT ~ 2 
keV and a 0.1-2.4 keV luminosity of 7.5 x 10 36 /i~ 2 W. 
From the PSPC image the emitting gas has a full extent 
of > 400 h^ 1 kpc, stretching north-east of the peak, 
which coincides with the position of NGC326. The best 
description of the radial profile they obtained is a combi- 
nation of a cluster /3-model and an unresolved component 
(possibly a hot galactic corona). They obtained a cluster 
core radius of 171 hr 1 kpc. Because of the presence of 



V../ YJ V IV^Ci'X VJ WallJlWll l Ll T | ^(J(J(J / 


00 h 58 m 22"6 


Ontiral nocntion ( Aionnn i 


4_26°51'58"3 


Redshift 


0.0477 


Distance 


141 ft -1 Mpc 


Magnitude (B T o) 


13.9 


1.4 GHz total flux density 


1.77 Jy 


1.4 GHz total radio luminosity 


10 24 6 ft" 2 W/Hz 


Integrated rotation measure 


-25 ± 5 rad m -2 


Overall spectral index 


0.8 


Radio source largest linear size 


180 h' 1 kpc 


Arcsec to kpc conversion factor 


l"= 0.63 h,- 1 kpc 



h = Ho/100 kms -i Mpc -i, 90 =0.5 
Wc use the convention S„ oc u~ a 



the compact hot corona centered on NGC326, Worrall 
et al. (1995) proposed that buoyancy effects of the radio 
plasma in the X-ray atmosphere could be responsible 
for the bending of the radio lobes. They suggested that, 
in projection, each radio lobe bends sideways into a 
tail as the backflow approaches the galactic corona. In 
order to explain the observed source shape, the buoyancy 
condition requires a backflow velocity for the lobes 
material of v ~ 3000 kms -1 . 

This is paper I of a series of papers dedicated to the 
radio galaxy NGC326. Here we give a description of the 
morphological, spectral and polarization characteristics 
of this radio galaxy. In Sec. 2 we summarize the obser- 
vations used in this paper. A description of the source 
morphology and its connections with the polarization 
and spectral properties is presented in Sec. 3. In Sec. 4 
and 5 we analyzed the jet collimation properties and the 
source physical parameters, respectively. A summary of 
the paper is given in Sec. 6. 

We will discuss the interpretation of the spectral and 
polarization properties respectively in the forthcoming pa- 
pers II and III. 



2. The VLA observations and data reduction 

A summary of the observations, including the VLA con- 
figuration, frequency, bandwidth, date and length of ob- 
servations, is reported in Table 2. 

The flux densities were brought to the scale of Baars et 
al. (1997) using 3C48 as primary flux density calibrator. 
The instrumental polarization and the polarization po- 
sition angle were calibrated using the sources 0116+319 
(0042+233 for 6 and 3.5 cm D array) and 3C138 respec- 
tively. The 3C138 polarization angle was assumed to have 
a value of -18° at L band and -24° at higher frequencies. 
Post calibration reduction was done using the National 
Radio Observatory (NRAO) AIPS package. The data for 
each observation were reduced using the standard self- 
calibration and cleaning procedures. The C array observa- 
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Table 2. VLA observations summary. 



Array 


Frequency 


Bandwidth 


Date 


Duration 




GHz 


MHz 




hours 


A 

A 


1.41/1.66 


12.5 


02-Dec-1984 


8.0 


A 


4.86 


12.5 


03-Dec-1984 


1.0 


B 


4.86 


25.0 


22-May-1985 


3.0 


B 


14.96 


25.0 


22-May-1985 


0.9 


C 


1.41/1.66 


50.0 


28-Sep-1985 


0.5 


C 


4.86 


50.0 


28-Sep-1985 


0.5 


C 


14.96 


50.0 


28-Sep-1985 


7.1 


D 


4.83 


50.0 


28-Sep-2000 


0.5 


D 


8.46 


50.0 


28-Sep-2000 


1.4 


D 


14.96 


50.0 


23-Nov-1985 


0.4 



Array 


Frequency 


Beam 






GHz 


arcsec 


mjy/beam 


A 


1.41 


l'/47 x 1724 


0.05 


A 


1.66 


1732 x 17 05 


0.05 


B 


4.86 


1719 x 1714 


0.029 


C 


1.41 


13782 x 127 77 


0.024 


C 


1.66 


11751 x 10771 


0.021 


c 


4.86 


3799 x 3777 


0.025 


c 


14.96 


1746 x 1725 


0.034 


D 


4.83 


14'.'08 x 12781 


0.06 


D 


8.46 


8797 x 7783 


0.05 


D 


14.96 


7731 x 4716 


0.08 



Col. 4: total intensity image noise. 



The angular size of NGC326 exceeds 4.4 arcmin. This 
represents a problem for high frequency imaging since the 
radio structure size is comparable to the field of view of a 
single VLA antenna. The old 2 cm observations suffered 
from this effect. Learning from experience, we followed a 
specific observation strategy with the more recent 3.5 cm 
observations. In order to reduce the problems due to the 
strong primary beam attenuation we splitted the obser- 
vations into two different pointings centered at appropri- 
ate positions on the wings of the source. The two data 
sets were reduced independently and the resulting couple 
of images was then combined using the task LTESS in 
AIPS. Table 3 summarize the relevant parameters of the 
single-array images. 

Data from different arrays were combined to improve 
uv-coverage and sensitivity. We combined the A and C ar- 
rays at 20 cm, the B, C and D at 6 cm, and the C and D 
at 2 cm. Each combined data set was self-calibrated. The 
images at different frequencies were cleaned and restored 
with the same beam using the AIPS task IMAGR. The 



"original" beams were very similar to the finally adopted 
ones. For the purposes of the spectral and polarization 
analysis we made data sets composed of matched resolu- 
tion images, using natural and uniform weighting. Finally 
we obtained four distinct resolution data sets with ex- 
actly the same cut in the maximum uv baseline at each 
frequency. Their resolutions are 10", 4", 2", and 172. 



Table 4. Combined array image parameters. The uv 
weighting function is indicated just after the frequency. 
UW and NW stand for uniform and natural weights, re- 
spectively. 



tion at 20 cm, already presented by Parma et al. (1991), 
has been reanalyzed for the purposes of this project. 

Table 3. Image parameters summary. 
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A+C 


1.41 uw 


0.22 


0.13 


A+C 


1.66 uw 


0.24 


0.14 


C+D 


4.86 nw 


0.07 


0.06 


D 


8.46 nw 


0.07 


0.01 


C+D 


14.96 nw 


0.20 


0.07 


4" X 4" beam 








A+C 


1.41 NW 


0.07 


0.04 


A+C 


1.66 nw 


0.08 


0.05 


B+C 


4.86 nw 


0.05 


0.03 


C+D 


14.96 nw 


0.04 


0.04 


2" X 2" beam 








A+C 


1.41 NW 


0.04 


0.03 


A+C 


1.66 nw 


0.05 


0.04 


B+C 


4.86 nw 


0.04 


0.03 


C+D 


14.96 nw 


0.03 


0.03 


1"2 X 172 beam 








A 


1.41 uw 


0.06 


0.06 


A 


1.66 uw 


0.07 


0.07 


B+C 


4.86 uw 


0.04 


0.05 


C+D 


14.96 uw 


0.04 


0.06 



Col. 3: total intensity image noise; 
Col. 4: polarization image noise. 



The equal resolution data sets are given in Table 4. All 
the images were corrected for primary beam attenuation. 
NGC326 is significantly polarized at all frequencies and 
resolutions. We produced images of the scalar fractional 
polarization (P v ), depolarization (DP^, defined as the 
ratio between the fractional polarization at v\ and i/ 2 ) 
and rotation measure (RM) . The polarization images were 
corrected for the non-Gaussian noise distribution of the 
polarized intensity. The RM was obtained using the AIPS 
task RM by a weighted fit of the position angle to the 
square wavelength at four frequencies. 

The spectral analysis has been performed with the pro- 
gram Synage++ (Murgia 2001). 
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Fig. 1. 1.4 GHz C array (contours) overlayed onto the optical image from the red Palomar Digitized Sky Survey. The 
contour levels are 0.5, 0.71, 1, 1.4, 2, 2.8, 4, 5.7, 8, 11, 16, 23, 32, 45, 64, 91, 130, 180 mJy (beam area)" 1 and the 
restoring beam is 13'.' 8 x 12'.' 8 with position angle —16.9° 

3. Source morphology, spectrum and polarization 

The observations described in the previous section allow 
us to investigate the morphology of NGC326 at increas- 
ing levels of angular resolution. In this section we describe 
the source morphology and its connections with the po- 
larization and spectral properties going from the extended 
structures (low resolution images) to the finest details 
(high resolution images). We refer to the different regions 
of the radio source as shown in Fig. 3. In particular, the 
two ellipses indicate the separation we marked between 
lobes and wings. 

3.1. Wings and plume 

The lower resolution images give us information on the 
extended structure of the source. At a resolution of about 
10" the lobes are slightly resolved and the most prominent 
components are the two wings, see Figs. 1 and 2. 



The wings bend and extend away from the lobe axis al- 
most by the same extent (2'). As already noted by Worrall 
et al. (1995), the overall 'Z'-shape symmetry of the source 
is broken by a low surface brightness plume located just 
above the end of east wing. This plume, evident only in 
the 20 cm and 6 cm images, does not follow the source 
symmetry. Fig. 1 shows a galaxy which is located, in pro- 
jection, in the middle of the plume. This corresponds to 
galaxy G6 of Werner et al. (1999) with a radial velocity 
difference of 1026 km s _1 with respect to the center of 
the cluster. If G6 is a cluster member it moves very fast 
relatively to the intra-cluster medium (slightly more than 
twice the line-of-sight velocity dispersion). Since the sound 
speed in the cluster is ~ 650 km s _1 , the motion of G6 
would be supersonic. Anyway, the association of G6 to the 
plume is unclear. 

Fig. 2 presents the 8.5 GHz map, with the E-ficld 
vectors superimposed (not corrected for RM, see below). 
The mean polarization percentage (10 "beam) in the 
wings is ~ 4% at 1.4 GHz and ~ 17% at 4.8 and 8.5 
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Fig. 2. 8.5 GHz D array image of NGC326. The contour levels are 0.15, 0.34, 0.69, 1.4, 2.9, 5.8, 12, 24 mJy (beam 
area) -1 and the restoring beam is 9" x 7'.' 8. The vectors lengths are proportional to the degree of polarization, with 
100 per cent corresponding to 20 arcsec on the sky, and their directions are those of the E-vector. The unresolved 
source at the tip of the east wing is a flat-spectrum background object. 
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Fig. 3. Source regions referred to in the text. 



GHz. Both wings show depolarization between 4.8 and 
1.4 GHz: the mean values of DP^g (10"beam) are 0.3 
and 0.2 for the east and west wing respectively. We cal- 
culated the RM at 10" resolution between the frequen- 
cies 1.4, 1.6, 4.8 and 8.5 GHz. We found a mean value 
of —25 ± 5 rad m~ 2 , which is consisten t with the galac- 
tic value ( Simard-Normandin et al. 1981 ), and a standard 
deviation ctrm = 52 rad m~ 2 . Given this value for the ro- 



tation measure, the apparent magnetic field direction at 
8.5 and 4.8 GHz is within few degrees from the real one. 

The apparent magnetic field is longitudinal and highly 
aligned along the main ridge of emission in both the east 
and, at least in the first part of, the west wing. At the 
wing edges the magnetic field bends to a circumferential 
configuration. 

The surface brightness of the plume is too low to give 
a reliable polarization measure. 

At a resolution of 10", the radio spectrum of the wings 
between 1.4 and 8.5 GHz can be computed out to their full 
extent. Fig. 4 shows a plots of the two-frequency spectral 
indices, a|;| and a|;|, as a function of the position along 
the two wings. The spectral indices have been sampled 
in circular boxes, of the same size as the beam, centered 
on the ridge lines of maximum emission. In the east wing 
there is a clear monotonically steepening of the radio spec- 
trum from the south lobe to the end of the wing: oq;\ and 
af'f increase from 0.6 and 0.7 up to 1.3 and 1.9, respec- 
tively. The west wing shows a different spectral behavior: 
c4;| and afj, increase respectively from 0.6 and 0.7, close 
to north lobe, up to 1.3 and 1.5, at a distance of ~ 60" 
(38 hr 1 kpc). In the remaining part of the wing, d\\ and 
a| | decrease and saturate to 1.0 and 1.3, respectively. 

The resolution of our images allows us to trace spec- 
tral index trends also along directions perpendicular to 
the wing ridge lines. There are significant lateral spectral 
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Fig. 4. Spectral index profiles along the wings ridge lines. The filled and open dots refer to and af f, respectively. 



index gradients with marked differences between the two 3.2. Lobes 
wings. 

The images at an angular resolution of 4" can be used to 

The detailed analysis of these spectral index profiles stud y the sha P e of the two lobes ( see Fi §' ^ Thc lobes 
will be presented in paper II are asymmetric in total emission, extent and dis- 

tance from the core (see also Fig. 6). The southern lobe 
has an ellipsoidal shape, while thc northern one is more 
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elongated and wider. The north lobe exhibits a shell-like 
feature which is particularly evident in the gradient im- 
age shown Fig. 6. This image, obtained with the AIPS task 
NINER, also reveals a sharp frontal border in the south- 
ern lobe and a straight edge in the east wing. By contrast 
the northern lobe is characterized by well defined lateral 
edges but has a smooth frontal border. Both lobes lack 
hotspots. The only structure with a significant brightness 
contrast with respect to the underlying lobe is the shell 
region in the north lobe. 

The mean polarization percentage at 8.5 GHz 
(10" beam) is 12% and 14% in the south and north lobe, 
respectively (see Fig. 2). Fig. 6 presents the 4.8 GHz im- 
age at 2" resolution, with the E-field vectors superimposed 
(not corrected for RM, see Sect. 3.1). The mean scalar 
fractional polarization at 4.8 GHz is 23% and 29% for 
the south and north lobe, respectively. The variation of 
the fractional point to point polarization at 4.8 GHz is 
strongly correlated with the enhancements seen in the to- 
tal intensity gradient image. In the lateral edges of the 
northern lobe and in the frontal head of the southern 



lobe, the fractional polarization is about 43%. In the 
sharp western edge of the east wing P4.8 reaches a value of 
- 50%. The mean values of DPH (10"beam) are 0.23 and 
0.17 in south and north lobe, respectively, while at 2" reso- 
lution they go up to 0.52 and 0.48. In both lobes, the mean 
RM at 10" resolution is about -20 rad mT 2 . As a conse- 
quence the Faraday rotation is negligible at 4.8 GHz and 
the apparent magnetic field in the lobes is circumferen- 
tial, consistently with the configuration deduced from the 
lower resolution 8.5 GHz map. The magnetic field config- 
uration is also circumferential around the shell structure, 
where the fractional polarization P4. 8 is about 30%. 

The mean fractional polarization polarization of the 
lobes as function of frequency and beam area is shown in 
Fig. 7. 

We have computed the average two-frequency spectral 
indices and aff g along slices perpendicular to the 
lobe axis from the lobe head back to the core using the 
2" resolution images. Slices are 3" wide, so that the spec- 
tral index measures are practically independent. The re- 
gions containing the core and the jet in the northern lobe 
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Fig. 6. Top panel: 4.8 GHz B+C array total intensity image. The contour levels are 0.14, 0.22, 0.36, 0.58, 0.93, 1.5, 
2.4, 3.9, 6.2, 10 mJy (beam area) -1 . The vectors lengths are proportional to the degree of polarization, with 100 per 
cent corresponding to 10 arcsec on the sky, and their directions are those of the E- vector. Bottom panel: 1.4 GHz A+C 
array gradient image (AIPS task NINER). In both images the restoring beam is 2" x 2". 
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South Lobe North Lobe East Wing West Wing Plume 



520cm ( m Jy) 


213 


576 


536 


288 


18 


8i x 62 (arcsec) 


30 x 22 


60 x 32 


100 x 50 


110 x 60 


90 x 60 


di x d'2 (kpc) 


22 x 14 


38 x 20 


63 x 32 


69 x 38 


57 x 38 




0.60 ± 0.04 


0.64 ± 0.02 


0.90 ±0.01 


1.04 ±0.01 


2.37 ±0.02 


(P1.4), 10"beam 


3% 


2.5% 


6% 


3% 




{Pi.s}, 10"beam 


13% 


15% 


22% 


16% 




(P 8 .5>, 10"beam 


12% 


14% 


17% 


14% 




(Pi.4>, 2"beam 


12% 


14% 








(Pi.s}, 2"beam 


23% 


29% 








(DPli), 10"beam 


0.23 


0.17 


0.27 


0.18 




{DPti}, 2"beam 


0.52 


0.48 








(RM) (rad m" 2 ) 


-21 ±6 


-18 ±4 


-34 ± 2 


-32 ±6 




(orm) (rad m~ 2 ) 


57 


52 


35 


65 





have been excluded from the statistics. Fig. 8 shows the 
spectral index trends in the lobes. In the south lobe a\'f 
is roughly constant around a value of 0.6 with a moderate 
increase to 0.7 at the lobe end; aff 9 increases from 0.7 to 
1.3. The spectral behavior of the north lobe is somewhat 
complex: ot\\ and aff 9 decrease respectively from 0.75 
and 1.35, at the lobe head, to 0.65 and 0.9, stay constant 
from 10" to 30", and then increase again up to 1.0 and 1.6 
in proximity of the core. The interpretation of the spec- 
tral profiles in the north lobe should be considered with 
caution since it is quite possible that the lobe, the jets and 
part of the west wing are overlapping each other because 
of projection effects. 

The main properties of radio lobes and wings are sum- 
marized in Tab. pi 



3.3. Jets 

Images of the jets at a resolution of 2" and 1"2 are 
shown in Fig. 6 and Fig. 9, respectively. The east jet 
can be traced out to a projected distance of 28.4/i _1 kpc 
until the south lobe, whereas the west jet extends out to 
7 kpc from the core, where it suddenly widens and 
merges into the lobe. Initially both jets are straight with 
a position angle p. a. ~ 43°. The east jet bends gradually 
(with a curvature radius of about 70") and approaches 
the south lobe almost horizontally. 

The fractional polarization P4.8 is about 14% in the 
east jet. The overlap with the lobe and the wing precludes 
estimation of the fractional polarization in the west jet. 
The apparent magnetic field is transverse in the east 
jet, whereas in the west jet it is initially transverse and 
becomes longitudinal in correspondence with the shell. 



a. 



30 



20 



o 

r-1 10 



(0 

a. 




Fig. 7. Mean fractional polarization in the lobes at dif- 
ferent angular resolution and frequency. The solid and 
dashed lines refer to the north and south lobe, respec- 
tively. 



We derived the spectral index between 1.4 and 4.8 GHz 
along the two jets as a function of the distance from the 
core using the 2" resolution images. The resulting trends 
are shown in Fig. 10. The spectral index stays almost con- 
stant at a value of 0.6 along both jets, although a signifi- 
cant scatter is present in the east jet. 

3.4. Radio cores 

The VLA images at maximum resolution combined with 
the optical Hubble Space Telescope (HST) image of 
NGC326 (Capetti et al. 2000) show that also the sec- 
ondary nucleus of the dumbbell hosts a radio core (see 
Fig. 9). According to the convention adopted by Werner 
et al. (1999) for the optical cores, hereafter we refer to the 
two radio cores as 'core 1' (the radio galaxy) and 'core 2' 
(the companion galaxy). We measure an angular distance 
between the two radio cores of 7"1 which corresponds to 
a projected separation of 4.8 h^ 1 kpc. The luminosity of 
core 1 is t ypical for a radio galaxy of this total power at 
1.4 GHz ( Feretti et al. 1984 ). Core 2 is an order of mag- 




Fig. 8. Spectral index profiles along the lobes at 2'.' resolution. The filled and open dots refer to a|;| and aff 9 , 
respectively. The shaded box indicates the region excluded from the spectral index calculation. 



nitude fainter and, at the sensitivity limits of our obser- 
vations, does not posses radio jets on kpc scale. 

The radio spectra of the two cores between 1.4 and 
14.9 GHz are shown in Fig. 11. Core 1 has a convex radio 
spectrum which peaks at about 10 GHz. Core 2 has a 
powerlaw radio spectrum with an index of 0.3. 



Table 6 summarizes the parameters of the two cores 
derived from the images at l'/2 resolution. 
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arcsec 



Fig. 9. 1.4 GHz A array (contours) overlaid with the HST image of the dumbbell galaxies (grayscale). The contour 
levels are 0.23, 0.33, 0.48, 0.71, 0.85, 1.1, 1.6, 2.3, 3.4, 4.8 mJy (beam area)" 1 and the restoring beam is T/2 x l'/2. 



Distance from Core (arcsec) 



Distance from Core (arcsec) 



Fig. 10. Profiles of the spectral index between 1.4 and 
4.8 GHz along the jets. 




v (GHz 



Fig. 11. Radio spectra of the two cores. Core 1 is the radio 
galaxy. 



4. Collimation and surface brightness of the jets 

We analyzed the evolution of jet collimation and surface 

brightness as a function of the distance from core. perpendicular to the jet axis for distances up to 40" and 

In order to quantify the variations of surface brightness 20" from the core, for the east and west jet, respectively, 
and width along the jets, we produced intensity profiles After having removed the baseline, the transverse bright- 
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Table 6. Core properties. 







core 1 


core 2 


1.4 GHz 


position (aj2ooo) 


00 ,l 58 m 22'.'6 


00' l 58 m 22'.'9 


1.4 GHz 


position (5j2ooo) 


+26°51'58'/7 


+26°51'52'.'7 


SlA GHz 


(mJy) 


5.33 ± 0.04 


1.16 ±0.04 


S4.8 GHz 


(mJy) 


10.42 ±0.03 


0.78 ± 0.03 


ail 




-0.5 


0.3 


PlA GHz 


(/i^WHz" 1 ) 


1.2 x 10 22 


2.7 x 10 21 


-P4.8 GHz 


(/i^WHz" 1 ) 


2.3 x 10 22 


1.8 x 10 21 



ness profiles are well fitted by a simple Gaussian over most 
of the jet length. 

We used both the 2" and l'/2 resolution images at 
20 cm. We measured the peak surface brightness 1$ and 
FWHM $0 and derived the deconvolved values I and $ 
by the first-order corrections: 



$ = ($n 



I = Jo ■ ($b/$ 2 + 1 



,1/2 



(Killcen, Bicknell & Ekers 1986) where $ b is the FWHM of 
the restoring beam. Fig. 12 displays the variation of $ and 
I as function of distance 9 along the jets. The steps of the 
transverse cuts are 1" and 0.6" in the 2" and 1'.'2 resolution 
images, respectively. The error bars represent the formal 
l-<7 confidence interval given by the fitting procedure and 
do not include the uncertainty introduced by the baseline 
subtraction. 

Apart for the initial region (0 <12"), where the jets 
are quite symmetric (see below) , the dependence of $ and 
I on are very different for the east and the west jet. 
The east jet spreads gradually up to 40" from the core, 
but has local regions that deviate significantly from the 
mean. Its deconvolved surface brightness drops by an or- 
der of magnitude in going from the core to a distance of 
15". Then, after a local knot of emission at =20", I 
decreases slowly for the remaining 20". On the contrary, 
the west jet expands abruptly at 12" from the core. Its 
deconvolved brightness initially decreases to a local mini- 
mum at =12" and then raises again peaking at =14". 
Beyond this peak the jet dccollimates originating the shell- 
like structure. 

The east/west width ratio within 12" from the core 
has a mean value of 0.98 and a range from 0.7 to 1.4. 
The surface brightness ratio has a mean value of 1.2 and 
a range from 0.49, at 2" 4 from the core, to 2, at 4"2 from 
the core. The overall jet symmetry for <12" suggests 
that the jet axes, at least at their base, lie in the plane of 
the sky. The asymmetry observed at larger distances may 
be due to strong projection effects and/or differences in 
the density of the gas surrounding the source. 

Both jets exhibit an average flat relation (I oc 

$ -1 - 35 , see Fig. 13). This behaviour indicates either that 
they depart significantly from the assumption of adiabatic 
flow, or they must be subject to a rapid longitudinal de- 
celeration (Fanti et al. 1982). For example, for a jet with 



no longitudinal component of magnetic field and a spectral 
index of 0.6 the adiabatic condition implies J oc <j>~ 3 - 4 Uj"" 3 , 
where Vj is the jet velocity; in this case the required de- 
celeration would be fj oc $-°- 68 . 



5. Physical parameters of the source 

Using the "minimum energy assumption" , we calculated 
with standard formulae (Pacholczyk 1970) the equiparti- 
tion internal energy (f/tot)j the energy density {u m - ln ) and 
the magnetic field (B eq ) for the jets and the extended 
components. We also calculated the equipartition pres- 
sure p C q = Umin/3. We assumed a volume filling factor of 
unity, equal energy in relativistic electrons and protons 
and a radio spectrum ranging from 10 MHz to 100 GHz 
with a spectral index of 0.6. 



5.1. Lobes and wings 

We assumed an ellipsoidal geometry for the lobes and a 
cylindrical geometry for the wings. Sizes for lobes and 
wings are those reported in Table 5. The global mini- 
mum energy parameters are listed in Table 7. We further 
calculated the variation of the equipartion parameters in 
ten regions in the wings (see Table 8 and Fig. 14). Worrall 
& Birkinshaw (2000) reported, for the cluster surrounding 
NGC326, a core radius of ~ 171 h^ 1 kpc and a pressure 
ranging from 7 x 10 -12 h 1 / 2 dyne cm -3 to 5 x 10~ 12 h 1 / 2 
dyne cm~ 3 at 100 kpc from the source. 



Table 7. Global equipartition parameters. 





Utot 


Umin 




Pcq 




Bcq 




erg 


erg/cm 3 




dyne/cm 


2 




N-Lobe 


6.2 x 10 56 


9.3 x 1(T 


12 


3.1 x 10" 


-12 


10.0 


S-Lobe 


3.4 x 10 56 


5.1 x 10" 


12 


1.7 x 10" 


-12 


7.4 


W-Wing 


3.3 x 10 57 


1.4 x 10" 


12 


4.8 x 10" 


-13 


3.9 


E-Wing 


3.1 x 10 57 


2.1 x 10" 


12 


6.9 x 10" 


-13 


4.7 



Therefore the inner radio lobes appear to be close to 
the pressure of the external gas while the wings appear to 
be under-pressured by a factor of 10. 

5.2. Jets 

We computed the minimum energy magnetic field and 
pressure in the jets using the deconvolved FWHM, $, and 
surface brightness, /, as derived in Sec. 4. We assumed for 
each segment of the jet a cylindrical shape with a cross 
section and length of 7T<1> 2 arcsec 2 and 1 arcsec, respec- 
tively. The minimum energy magnetic field varies from 40 
to 9 [iG going from the core outwards. The equipartition 
pressure varies from 5 x 10 -11 to 3 x 10~ 12 h 1 / 2 dyne/cm 2 
(see Fig. 15). 

From the X-ray luminosity reported by Worrall et al. 
(1995) for the hot galactic corona (i.e. 10 41 hr 2 erg/sec), 
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Fig. 12. Top: deconvolved jets FWHM $ (left panel) and brightness / (right panel) profiles from the 2" resolution 
20 cm image. Bottom: deconvolved FWHM (left panel) and brightness (right panel) ratios from the l'/2 resolution 
20 cm image. 



we computed the expected external pressure on the jets. 
We assumed for the corona a temperature of 1 keV and a 
gas density profile n = n - (l + (r/r c ) 2 ) , where n is the 
central density while r and r c are the distance from the 
galaxy and the corona core radius, respectively. We found 
that the jets are roughly in pressure equilibrium with the 
surrounding gas (see Fig. 15) for a corona core radius r c = 
3 hr x kpc and b = 3; in this case the central density results 
no ~ 0.2 h 1 / 2 cm~ 3 . A model with r c = 2 h^ 1 kpc, b = 2 
and no — 0.3 h 1 ^ 2 cm -3 gives an equivalent result. 



6. Summary 

We have presented multi-frequency data of the radio 
galaxy NGC326, together with previously published data 
that were re-analyzed for the purpose of this work. These 
data allowed us to investigate the source morphology at 
different levels of spatial resolutions. 

The source structure is complex at all resolutions and 
may be dominated by projection effects. 

At the lowest resolution the twin wings appear fairly 
symmetric but the 'Z '-shape is disturbed by a plume of low 
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Table 8. Wings cquipartition parameters (a = 0.6). 

distance u min P cq B oq 

arcsec erg/cm 3 dyne/cm 2 fiG 
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Fig. 13. Deconvolved brigthness as a function of the de- 
convolved FWHM, /(<£>), for the west (open dots) and east 
(filled dots) jet. The line has a slope of -1.35. 



surface brightness and very steep spectral index (a > 2). 
The plume overlap, in projection, a possible fast-moving 
member of the galaxy cluster. 

At the intermediate resolution of about 4" the source 
lobes appear asymmetric in shape, size, brightness and 



projected distance from the core. No compact components 
(hostspots) are present in the lobes. 

At higher resolution, the jets are very symmetric in 
both flux density and width within the inner 12" from 
the core. This indicates that the axis of the jets, at least 
at their base, lies on the plane of the sky. Beyond this 
distance, the east jet bends gradually with a curvature 
radius of about 70" reaching the south lobe. The west 
jet is straight but, after recollimation expands abruptly, 
giving rise to a shell-like structure embedded in the north 
lobe. 

Overlapping our highest resolution VLA images of the 
core region with the HST image of the dumbbell, we found 
that also the companion of the radio galaxy hosts a radio 
core. 

NGC326 is significantly polarized at all frequencies 
and resolutions. The magnetic field configuration is 
circumferential in the lobes and highly aligned with the 
ridge of emission in the wings. 

We traced the profile of the spectral index between 1.4 
and 8.5 GHz for the entire wing length. The spectral index 
increases gradually going from the regions in proximity of 
the lobes to the wing ends. The spectral index distribution 
in the lobes is complicated, presumably due to the mixing 
with the wings. Both jets show a constant spectral index 
profile with a mean value of 0.6. 

Finally, we computed the minimum energy, magnetic 
field and pressure in wings, lobes and jets. We found that 
the lobes are close to equilibrium with the external gas 
pressure while the wings are under-pressured by a factor 
of 10. 

These observations allow us to investigate in detail the 
spectral and polarizarion properties of this peculiar object 
at different levels of spatial resolution and constitute an 
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Fig. 14. Minimum energy magnetic field (solid line, left y- 
axis) and pressure (dashed line, right y-axis) in the wings. 
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Fig. 15. Minimum energy magnetic field (solid line, left 
y-axis) and pressure (dashed line, right y-axis) in the jets. 
The dotted line represents the expected thermal pressure 
from the cluster gas and from the hot galactic corona 
model described in the text. 

excellent data set which permits testing the "standard" 
theory of the synchrotron aging analysis. 

The interpretation of the data will be discussed in pa- 
pers II and III. 
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